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We study some implications of the presence of two new scalar weak doublets beyond the stan¬ 
dard model which have zero vacuum expectation values and are charged under an extra Abelian 
gauge symmetry. The additional gauge sector does not couple directly to standard-model parti¬ 
cles. We investigate specifically the effects of the scalars on oblique electroweak parameters and on 
the interactions of the 125 GeV Higgs boson, especially its decay modes h —> 77,7 Z and trilinear 
self-coupling, all of which will be probed with improved precision in future Higgs measurements. 

Moreover, we explore how the new scalars may give rise to strongly first-order electroweak phase 
transition and also show its correlation with sizable modifications to the Higgs trilinear self-coupling. 


I. INTRODUCTION 


The recent discovery |l|, l2[ at the Large Hadron Collider (LHC) of a Higgs boson with mass around 125 GeV and 
other properties consistent with the expectations of the standard model (SM) serves as yet another confirmation that 
it is a remarkably successful theory. Nevertheless, it is widely believed that new physics beyond it is still necessary at 
least to account for the compelling experimental evidence for neutrino mass and the astronomical indications of dark 
matter Q. 

Among a great many possibilities beyond the SM are those with enlarged scalar sectors. Scenarios incorporating 
a second Higgs doublet are of course highly popular in the literature 0, [5]. Of late models with three scalar weak 


doublets have also been gaining interest [6 14|, as they can provide dark matter (DM) candidates 13( and/or an 
important ingredient for the mechanism that generates neutrino mass 14]. 

Here we consider this three-scalar-doublet possibility, particularly that in which two of the doublets possess zero 
vacuum expectation values (VEVs). The theory also involves a new Abelian gauge symmetry under which these two 
doublets are charged, while SM particles are not. As a consequence, the extra scalar particles do not couple directly 
to a pair of exclusively SM fermions. Because of the absence of their VEVs and couplings to SM fermion pairs, these 
scalars have been termed inert in the literature Q. However, being members of weak doublets, these scalars have 
interactions with SM gauge bosons at tree level. In addition, the gauge boson associated with the new gauge group is 
taken to have vanishing kinetic mixing with the hypercharge gauge boson. Accordingly, the additional gauge sector 
can be regarded as dark. 

With these choices,the scalar sector of the theory corresponds to one of the three-scalar-doublet models catalogued 
and studied in Ref. [8( in terms of all possible allowed symmetries. In the present paper, we entertain the scenario 
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described above and explore some implications of the presence of the inert scalars. Specifically, we analyze constraints 
on them from collider measurements on the Higgs boson and from electroweak precision data. In addition, we look at 
the potential impact of the scalars on the Higgs trilinear self-coupling, anticipating future experiments that will probe 
it sufficiently well. To evaluate the coupling, we will employ the Higgs effective potential derived at the one-loop level. 
Moreover, we examine how the new particles, which we choose to have sub-TeV masses, may give rise to strongly 
first-order electroweak phase transition (EWPT), which is needed for electroweak baryogenesis to explain the baryon 
asymmetry of the Universe. As it has been pointed out in the context of other models that the strength of EWPT 
could be correlated with sizable modifications to the Higgs trilinear self-coupling 15our results will indicate how 
this may be realized in the presence of the new doublets. 

Due to their tree-level interactions with SM gauge and Higgs bosons, the lightest of the inert scalars cannot serve 
as good candidates for DM, as they annihilate into SM particles too fast and hence cannot produce enough relic 
abundance. To account for DM, one needs to have a more complete theory, but we assume that the additional 
ingredients responsible for explaining DM have negligible or no impact on our scalar sector of interest, so that they 
do not affect the results of this paper . 1 

The plan of the paper is as follows. In the next section, we describe the scalar Lagrangian and address some 
theoretical constraints on its parameters, especially from the requirements on vacuum stability. Since the extra scalar 
doublets couple to the standard Higgs and gauge bosons and include electrically-charged members, they contribute 
at the one-loop level to the Higgs decays h —> 77 and h —> 7 Z which have been under intense investigation at the 
LHC, the former channel having also been observed. We determine their rates in section mu where we also start 
our numerical analysis by exploring the charged scalars’ impact on these processes. In section Hvl we calculate the 
contributions of the new doublets to the oblique electroweak observables S and T, on which experimental information 
is available. Sections IVl and fVTl contain our treatment of the new scalars’ effects on the trilinear Higgs couplings and 
on the electroweak phase transition, respectively. After deriving the relevant formulas, we perform further numerical 
work in these sections. In section lVlII we discuss additional results and make our conclusions after combining different 
relevant constraints. A few appendices contain more discussions and formulas. 


II. SCALAR SECTOR 


A. Lagrangian 


Compared to the SM with the Higgs doublet <f>, the scalar sector is expanded with the addition of two weak doublets, 
774 and 772 . The theory also possesses an extra Abelian gauge symmetry, U(1 )d, under which 7^2 carry charges +1 and 
— 1, respectively, whereas SM particles are not charged. Accordingly, one can express the renormalizable Lagrangian 
for the interactions of the scalars with each other and with the standard SU(2)i x U(l)y gauge bosons, 1 ^ 2,3 and 
B , as well as the U(1 )d gauge boson C, as 

£ = + (V^V^rn + (V^g 2)^72 - v, (1) 


where the covariant derivative V^ = cU + ( ig/2)rjW+ igyQyB^ + igoQcC M also contains the gauge couplings g , 
gy 1 and <?d, Pauli matrices tt, 2 , 3 , and U(l)y i £> charge operators Qy,c, while the scalar potential is 


V = ^$ 1 $ 


1A1 vim 


1 A 2 V 2 V 2 


+ + \\2\{g\gi) 2 + 


< \ 3 i$' i '<I>7[77i + A 32 < h^ < I> 77^772 + A4i<I> t 77i77j + A42$ t 77277|<h 

A 5 $ t ? 7 i $ t 772 + + Ag 77^774 772 772 + A 7 t 7 + 77277+ 774. 


( 2 ) 


1 In a recently proposed scotogenic model El , the inert scalars participate in the mechanism to generate light neutrino masses via one-loop 
interactions with new fermions which include good DM candidates. In such a context, our results would likely be modified. 
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Thus Qc$ = 0 and Qcih{v 2 ) = +Vi(~ r l 2 )- The parameters /r 2 2a an d ^i, 2 a, 3 o, 4 a, 6,7 with a = 1,2 are necessarily 
real because of the hermiticity of V, whereas A 5 can be rendered real using the relative phase between $ and 
Assuming that the U(1 )d symmetry stays intact, after electroweak symmetry breaking we can write 


* = ( ^ (T ,° + h) ) ’ ’ 7o= (^°)’ V2r, 0 a =Rer, 0 tt +ilmT, 0 a , (3) 

where h represents the physical Higgs boson, v ~ 246 GeV is the vacuum expectation value (VEV) of $, and H+ and 
77 ° denote, respectively, the electrically charged and neutral components of ry a , which has no VEV. 

From the terms in V that are quadratic in the fields, it is straightforward to extract the mass eigenstates of the 
scalars. Thus the masses of h and H± 2 at tree level are given by 

™h = l* ? + |AiU 2 , m 2 H ± = ula + \\ 3 aV 2 . (4) 

The A 5 part in eq. ED causes mixing between the electrically neutral components 77 ° and 77 °*, which are then related 
to the mass eigenstates xi and X 2 according to 



tan(20) 



_ X5V 2 _ 

2m^ ± - 2 m 2 H± + (A 42 - A 41 ) v 2 


sg = sin#, 


the resulting eigenmasses being given by 


\ (m 2 H ± + m 2 H ±^ + \ (A 4 i + A 42 )' 


T 


2 „± - m 2 H± + \ (A 42 - A 4 i) v 2 


H. 


+ |A§^ 


(5) 


( 6 ) 


Hence the U(1 )d charges of xi ,2 are the same as (opposite in sign to) that of 771 (r) 2 ) and m Xl < m X2 . 
Alternatively, instead of Xa, one can choose to deal with their real and imaginary parts, 


S a = V 2 Rexo, V a = V 2 Imx'a, 


(7) 


which are CP-even and CP-odd states, respectively, and share mass, ms a = rn-p a = 77i Xa . From eq. ED, one then has 
in matrix form 


^ Re77i ^ 


/ 

eg 

sg 

0 

° ^ 


( s 1 \ 

Re 77 ° 



-sg 

eg 

0 

0 


s 2 

Im77i 



0 

0 

eg 

sg 


Pi 

V Im7 ?2 ) 


l 

0 

0 

sg 

-eg y 


\v 2 ) 


( 8 ) 


where the mixing matrix is orthogonal. 

Later on, to simplify the analysis, we will concentrate on the scenario in which A 5 is negligible compared to 
the other A’s in V. In that case, as eq. ED indicates, the 771 - 772 * mixing is small, 9 <C 1, provided that A 5 V 2 -C 
2m 2 ± — 2m 2 ± + (A 4 i — A 42 ) v 2 . Furthermore, one can see from eq. ED that at the same time xi and X 2 can be close 

Hi H 2 

in mass if hXsv 2 -C 


m 2 „ ± - m 2 H± + \ (A 4 i - A 42 ) v 2 




<C m : 


Xl ‘ 


B. Theoretical Constraints 


The parameters of the scalar potential are subject to a number of theoretical constraints. The stability of the 
vacuum implies that V must be bounded from below. As shown in appendix [Aj with A 5 being negligible, this entails 
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that for a = 1, 2 


+ 


Ai > 0, A2 a > 0, A3 a + A° 0 + \J A1A20 >0, Ag T A? + VA21A22 > 0, 

\/ A1A2I A22 + \/Al (Aq + Ay) + \/ A2I (A32 + A42) + \/ A22 (A31 + A4-J 

,1 !/2 


2 (\A^ 


+ A, 


31 


+ A° 


41 


) (v^iAa 


+ A; 


32 


+ A° 


42 


A 21 A 2 


+ Ag + Ay ^ 


> 0 , 


(9) 


where A° = Min(0, A^). 

The /i 2 and A parameters in V also need to have such values that its minimum with the VEV of $ {rj a ) being 
nonzero (zero) is global. This is already guaranteed j<J] by the positivity of the mass eigenvalues in eqs. 0 and 0. 

In addition, the perturbativity of the theory implies that the magnitudes of the A parameters need to be capped. 
Thus, in numerical work our choices for their ranges, to be specified later on, will meet the general requirement 
| A* | < 87 t, in analogy to that in the two-Higgs-doublet case [20]. 


III. RESTRICTIONS FROM COLLIDER DATA 


The kinetic portion of the Lagrangian in eq. 0 contains the interactions of the new scalars with the photon and 
weak bosons. 


C D tH+o' J H- (eA„ - g L Z„) + H+H- (eA - g L Z) 2 


+ 


W 

2c w 


(x*< 


C28 ( Xld^Xl ~ X*2d >J 'X2 


) + s 2 e (xt 9^X2 + X&Xl) 


ry .J ry Z 

4 C 2 XaXaZ 


ce 


(H+d^xX + H+d» X 2) + s e [H+d^X 2 - H+d» X 1 ) 


+ f { L 

+ {H+H, + xlXa) W+»W~ , 


W-- H.C.} 


where summation over a = 1,2 is implicit, 

Xd^Y = Xd+Y - Yd^X 


9L = ^(2sl-l), 


( 10 ) 


( 11 ) 


c w = cos0 w = (1 — s^,) 1 / 2 , with 0 W being the usual Weinberg angle, C 20 = cos(20), and S 20 = sin(20). One can 
alternatively write eq. m in terms of the real and imaginary components S a and V a of \a , which becomes more 
lengthy and is relegated to appendix [B] 

We now see that data from past colliders can lead to some constraints on the masses of the new scalars. Based on 
eq. cdd. we may infer from the experimental widths of the W and Z bosons and the absence so far of evidence for 
nonstandard particles in their decay modes that for a, b = 1,2 


m H ± + m Xb > mw , 2m H ± > mz, m xa Ym Xb > mz■ (12) 

The null results of direct searches for new particles at e + e _ colliders also imply lower limits on these masses, especially 
those of the charged scalars. 2 For these reasons, in our numerical work we will generally consider the mass regions 
m Xa — 50 GeV and > 100 GeV. 

In addition to the requirements in the preceding paragraph and the vacuum stability conditions in eq. ©, when 
selecting the inert scalars’ parameters we take into account also the Higgs mass which will be estimated at the 


2 A recent investigation [ 2 H concerning the effects of the corresponding particles in the simplest scotogenic model [531 on the relevant 
processes measured at LEP II suggests that such charged scalars may face significant constraints if their masses are below 100 GeV. 
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one-loop level in section 0 and then limited to = (125.1 ± 0.1) GeV, well within the ranges of the newest 
measurements 


23 


24]. More specifically, we will therefore make the parameter choices 


0 < A 2a , |A 3q | , |A 4 „|, |A 6 |, |A 7 | < 3, j/z| a | < (800 GeV ) 2 
IAg| < O.OlMin(A 2a , |A 3 a |, |A 4 a |, |Ag |, |A 7 |). 


(13) 


The recently discovered Higgs boson may offer a window into physics beyond the SM. The presence of new particles 
can give rise to modifications to the standard decay modes of the Higgs and/or cause it to undergo exotic decays Ehl . 
As data from the LHC will continue to accumulate with improving precision, they may uncover clues of new physics 
in the Higgs couplings or, otherwise, yield growing constraints on various models. Here we address some of the 
potential implications for our scenario of interest. Especially, the existing experimental information on the possible 
Higgs decay into invisible/nonstandard final states 26 -3c| and on the observed h —> 77 mode 24, 31 1 can supply 
further restrictions on the inert scalars. 

The Higgs boson couples to a pair of them according to 

£ D ^ [(/4 - m %J H+ iff + (/4 - m 2 H2 ) H+ 

+ { 4 t 4 1 + 4 n L - "* 2 i) x*xi + {4^22 + 4^21 - m\ 2 ) V2X2 

+ cgs e (/x21 - m 22 ) (X1X2 + X2X1)] , (14) 

from the V part of eq. ©• In view of the mass choices made above, it follows that the decay modes h —> XaXb: if 
kinematically allowed, contribute at tree level to the total width of the Higgs boson and are the leading channels into 
nonstandard final states in the model. Their rates have the form 




ml — rn ‘\a. _ m Xb ) _ '^ m Xa m Xb'> 


(15) 


where 


2 2 1 2 2 2 

^XiXi — C 6^21 s el 1 22 — m Xl , 

Cxlx 2 = C x ix\ = °s s s (/x 2 i — M22) ■ 


rx 2 2,2 2 ^2 

Cxix 2 — c e ^22 + s e^ 2 i ~ m X2 , 


The combined branching ratio of these decays is 

B(^X*X') = 


J 2 a,b T ( h ^ X*aXb) 
r h M + Ea,6 r ( h “»■ XlXb) ’ 


(16) 


(17) 


where T® M is the SM Higgs total width and only channels satisfying m Xa + m Xb < rrih contribute to the sums. 
Numerically, we adopt T® M = 4.08 MeV (32| corresponding to = 125.1 GeV. If these channels are open, we will 
require B(h —► XlXb) < 0.19, based on the latest analysis of the Higgs data 


The potential impact of the inert scalars can also be realized through loop diagrams. Of much interest are their 
contributions to the standard decay channels h 77 and h —> 7 Z , which are already under investigation at the 
LHC. In the SM, they arise mainly from top-quark- and VE-boson-loop diagrams. These modes receive additional 
contributions arising from the diagrams drawn in figure [TJ with vertices from eqs. (HOD and (1141) . 3 Their 


3 At the one-loop level, the charged (charged and neutral) inert scalars also induce h —»• 7 C ( h —>■ ZC , CC ) involving the massless dark 
gauge boson C. These decay modes may be challenging to detect with C being invisible, as their rates are expected to be roughly of 
similar order to those of the 77 and 7 Z channels. 
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decay rates are readily obtainable from those in the case of only one inert doublet 33|. Thus we get 

2 


T(h 77 ) = 


T(h -A 7 Z) = 


a 2 GF , m'l 

128v / 2tt 3 


7 ^ 1/2 ( K *) + ^i 7 ( k vk) + 


2 m 2 _ „2 

m ff ± 7*2(3 


a 77 




( K ffa±) 


(18) 


aG^rri^ (m 2 h - m 2 z Y 


647t 4 to 3 


6 - 16s?„ 


3c„ 


^1/2 Ct) + c wA 7 (kwaCw) 


2 _Lm 2 ff± -,& 


1 ~~ 2s w y- "" H „ 

C w . nzlb 

a=l 


^0 ( K £r±>CiJ±) 


(19) 


where a = <? 2 s^/( 47 r) is the fine-structure constant, the expressions for the form factors A 77 ’/ 7 ^ are available from 
ref. 34 j, the Aq 7,7Z terms originate exclusively from the Hf 2 diagrams, kx = 4771^/m 2 , and (x = 4 m 2 x /m 2 z . 


^wv 7 


■<: h u 


' -ww 7 > ^ 




FIG. 1: Feynman diagrams for the contributions of the new charged scalars Hf 2 to the Higgs boson decays h 
h —y 7 Z. The triangle diagram with the gauge boson legs interchanged is not shown. 


77 , and 


We can already test the new contributions to h —y 77 , which has been observed at the LHC, unlike the 7 Z channel. 
For the 77 signal strengths, the ATLAS and CMS Collaborations measured ct/osm = 1-17 ± 0.27 [iH and 1.13 ± 0.24 
24j . respectively. These numbers need to be respected by the ratio of T(/i —y 77 ) to its SM value, 

T(h -y 77 ) 


~ T(h 


77) 


SM 


Its 7 Z counterpart, 


= 


T{h -y 7 Z) 
T(h -y ryZ) SM 


( 20 ) 


( 21 ) 


will be probed by future experiments. 

To illustrate the effects of the inert scalars on 7 ?. 77 and IZ-yz, and possible (anti)correlation between them, we display 
in figure [5] the distribution of 5000 benchmark points on the (TZ^z ^TZ-yy) plane which satisfy the vacuum stability 
requirements in eq. m, the constraints from W and Z decays in eq. m, and the parameter limitations in eq. m 
We notice that many of the 1Z 11 values are close to 1 and within the allowed ranges from ATLAS and CMS. The 
plot also reveals that for the 1Z 71 points compatible with the LHC data the values of T(h —y 7 Z) do not differ from 
its SM value by more than 10% or so. Furthermore, there is a positive correlation between 7 £ 77 and 7 Z 7 z, which 
is much like the situations in a different recent model with two inert doublets [gj and in the case of only one inert 
doublet |33l l35l 1361 ]. This can be checked experimentally when the 7 Z mode is observed in the future. 


IV. ELECTROWEAK PRECISION TESTS 


The interactions of the new doublets with the SM gauge bosons described by eq. m bring about modifications, 
AS and AT, to the so-called oblique electroweak parameters S and T which encode the effects of new physics not 
directly coupled to SM fermions 371. At the one-loop level [3, 37] 


a AS Azz (mz) ~ Azz(0) 


• s 2 

^W°W 


4 c 2 
aAT = 


- A' ( 0 ) - 


c w s w 


A 7 z( 0 )> 


Aww(0) Azz(0) 


w 


m% 


( 22 ) 
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FIG. 2: The effects of the new charged scalars on the ratios of the rates of Higgs decay channels h 
to their respective SM values for 5000 benchmark points as described in the text. The blue point marks 
region between the green (magenta) horizontal lines represents the one-sigma range of the ATLAS (CMS) 


—>• 77 and h —► 7 Z 
the SM value. The 
data @, 0 . 


where the functions Axy (q 2 ) can be extracted from the vacuum polarization tensors A^ Y (q 2 ) = A xy (q 2 ) g' u ' + 
[q^q v terms] of the SM gauge bosons due to the new scalars’ loop contributions, and A' XY ( 0) = [dAxY (q 2 ) /dq 2 ] q 2 =0 . 
In our numerical analysis below, we will impose 


AS = 0.05 ±0.11, AT = 0.09 ±0.13, (23) 

which are based on the results of a recent fit 0] to electroweak precision data for a Higgs mass rrih = 125 GeV. 

The contributions of the inert scalars to AS and AT arise from the diagrams depicted in figure [3] After evaluating 
them, we arrive at 4 


where 


as = s 


ln + , ~ 5ra k * 5m h 


TTL iri TTl rr A 


6 (m^ - m| 2 ) 

( m xi + m 2 X2 ) (m 4 Xi - 4 m^m 2 , + m 4 J 


2 v-x 

s 29 


(ml, — m 


XI 


2 y 

X 2 / 


ln ■ 


‘Xi 


m y 


AT = 


1 


8an 2 v 2 


c 2 s T (m H ± ,m Xl ^+ c 2 0 T (m H ±, m X2 ) + s 2 0 T (m H ±, m x 
+s 2 g T ( 'm H ±,m Xl ) - 4cgSgT(m Xl ,m X2 ) 


T(m, n) = 


m 2 ± n 2 


2 2 2 

m n m 

-111-. 

2 2 2 

m z — n z 


(24) 


(25) 


(26) 


In figure |4] we present the distribution on the (AS, AT) plane of the inert scalars’ contributions for the 5000 
benchmarks employed previously for figure [2] Evidently, it is possible for the masses of the charged scalars to be as 
small as 100 GeV and still be compatible with the electroweak precision measurements. However, we find that the 
lighter one of the inert neutral scalars, xi> must be heavier than about 90 GeV, which is a stronger condition than 
that inferred from the LEP constraint on the invisible width of the Z boson. This also makes the bound from the 
data on the Higgs invisible/nonstandard decay irrelevant. 


4 Their counterparts in the case of only one inert scalar doublet were computed in Ref. l39fl ■ 



















7 vA/Vl 
\ 


kA/v s 7? Z 7vVWVv*7j^ Zwv^ \as^Z 

f B+ H+ 2 

W+ aa /1 i/vv^VF -1 " VF^ vWi j/vV'T'F -1 ' 

Ya X'2'l 


ZvW J/W*^ ^A^VWZ 
\ / 

X 24 
Ha, Xa 

x I 

W+^V<7<A^W + 


FIG. 3: Feynman diagrams for the contributions of the inert scalar doublets to the oblique electroweak parameters AS and 
AT. 



FIG. 4: The contributions of the inert scalar doublets to the oblique electroweak parameters AS and AT for the 5000 benchmarks 
used previously. On the left panel, the palette belongs to the lighter neutral inert scalar’s mass, m xi , in GeV. On the right 
panel, the palette belongs to the lighter charged scalar mass, mH x , in GeV. The different contours represent 68%, 95%, and 
99% confidence level, respectively. The blue point at (0,0) marks the SM value. 


V. HIGGS TRILINEAR COUPLING 


Since the new scalars couple directly to the Higgs boson, their presence can cause its trilinear self-coupling, A hhh, 
to shift from its SM prediction. Such a modification could translate into detectable collider signatures, especially at 
a future e + e _ machine such as the International Linear Collider [401 ] where the coupling can be measured with 20% 
precision or better at a center-of-mass energy yfs = 500 GeV if the integrated luminosity is 500 fb -1 . 

To derive the formula for the mass-dimension Higgs trilinear self-coupling in the presence of extra heavy particles, 
we follow the steps taken in ref. [41|. It is just the third derivative of the Higgs effective potential, namely 


^hhh — 


d 3 


Kff = °(7> 


(27) 


where <p is the classical Higgs field and u ( ( 7 is the potential evaluated at temperature T = 0. We estimate the 
potential at the one-loop level in the so-called DR' scheme [ 42 I, 0| where it has the form 


U T=0 

V eS 


/ \ £*12 
m = 


Al 4 

■ 


E 


Tli 


(™?(7 >) 2 ( ln m2 i(v) 


647r 2 


A 2 


(28) 


In the sum above, the index i runs over all the contributing particles, nt stands for the number of internal degrees 
of freedom of the ith particle, with a minus sign added if it is a fermion, mf((p) is its field-dependent squared mass, 
and A is the renormalization scale which we choose to be the Higgs mass, A = 125.1 GeV. More explicitly, nh = 1, 
ng = riz = n 7 = 3, n\y = 6, n t = —12, and n Xa = n H ± = 2, where Q refers to the Goldstone bosons. We have 
collected the formulas for the various relevant m 2 (tp) in appendix [Cl 
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At tree level we have p\ = — X\v 2 /2 = jl\, but it receives the one-loop correction 


nimfmf 




(29) 


which follows from dV^~°(p)/dp = 0 set at p = v ~ 246 GeV, where mf = m?(p) and mf = dm?/dtp. Then the 
Higgs mass at the one-loop level, which is nothing but the second derivative of is given by 


ml = Xiv 2 + 


32t r 2 


mf mf - 


+ (m 2 ) 2 ) In— m?ml + 


p—v 


(30) 


where the first term is the familiar tree-level contribution, the second term is the radiative one-loop correction, and 
to? = d 2 m 2 / dp 2 . Accordingly, with mf being fixed to its empirical value, as Ai is varied along with the other scalar 
couplings it can be bigger or smaller than its tree-level value Ai = mf/v 2 ~ 0.258, depending on the size and sign of 
the loop contribution in eq (l30l) . 

Incorporating eq. (I3(J1) into eq. (1271) . one then obtains 


Xhhh 


3 m\ 


V 





•m2 2 

m- mf 




tp—V 




(31) 


where mf = d 3 m? / dp 3 . Its SM counterpart, Afff, has the same formula, except that in the sum i runs over SM fields 
only. 

According to eq. (13T1) and appendix [€] the Higgs trilinear self-coupling is a function of the couplings A 3 a + A 4 a and 
A 3 a of the inert neutral and charged scalars, respectively, to the SM Higgs doublet, i.e. through the field-dependent 
masses and their derivatives. Since A 3 ai 4 a are related to the scalars’ physical masses via eqs. © and ©, the Higgs 
trilinear coupling also depends on them. To illustrate how the inert scalars’ couplings and masses affect A hhh, we 
define the relative change 


A = 


A hhh A S ,M 


\SM 

A hhh 


(32) 


with respect to the SM prediction. Then in figure [5] we graph A versus | A 32 + A 42 1 and | A 32 1, respectively, for the 5000 
benchmark points employed earlier. On the same plots we also show the mass distributions of the inert neutral and 
charged scalars, respectively. 

It is clear that in the presence of the inert doublets the trilinear Higgs self-coupling can be enhanced or reduced by 
up to roughly 150% relative to the SM contribution to it. One realizes that, for either large or small (charged and/or 
neutral) scalar masses and couplings to the SM Higgs doublet, this enhancement or reduction of the trilinear coupling 
is the effect of the superposition of different contributions which could be constructive or destructive. 

The new scalars’ impact can be further seen in figure El which illustrates their loop effects. Specifically, it displays 
the relative changes of the trilinear Higgs coupling, the Higgs mass, and the parameter due to radiative corrections 
versus the Higgs quartic self-coupling Ai, where 


^A hhh — A hhh 3Al77, 8m^ — TO^ AiV , 


(33) 


5p\ is defined in eq. (l29l) . and = (i\ + 5p\. 

We remark that the Higgs quartic self-coupling, which at tree level is defined by the Higgs mass, can have a wide 
range from about 10 -4 to 0.5. This is due to the fact that much of the Higgs mass arises radiatively, as the right 
plot in figure El indicates. More precisely, to^ can be fully radiative for small Ai values or get a negative radiative 
correction for large Ai values, those greater than its tree-level one, Ai. One can see from the top-left and bottom plots 
in the figure that similar remarks could be made concerning Xhhh and In particular, each of these parameters 
may be fully radiative for small Ai and also can receive radiative corrections which are negative. 
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FIG. 5: The changes of the Higgs trilinear self-coupling relative to its SM value versus the absolute values of the SM Higgs 
doublet couplings, A 32 + A 42 and A 32 , to the heavy neutral (left) and the heavy charged (right) scalars, respectively. On the 
palettes, we read the heavy neutral (left) and charged (right) scalar masses in GeV. 


VI. ELECTROWEAK PHASE TRANSITION 


It is well-known that one of the reasons why the SM fails to produce successful baryogenesis 44[ is the fact that the 
EWPT is not strong and consequently cannot suppress processes that violate the conservation of baryon plus lepton 
numbers, B+L , in the broken phase [4a]. The suppression of anomalous B+L -violating processes in the broken phase 
happens if the criterion for strongly first-order EWPT 00, 

v c /T c > 1, (34) 


is fulfilled, where v c is the Higgs VEV at the critical temperature T c at which the effective potential exhibits two 
degenerate minima, one at zero and the other at v c . Both T c and v c are determined using the full thermal effective 
potential 48, 0 


Veg((p,T) = V3f°(<p) 


(35) 


at a finite temperature T, where 


Jb.f (r) = dx x 2 In 


1 =F exp - 


\/x 2 + rj 


(36) 


the upper (lower) sign referring to a boson (fermion). To V e ff(tp,T) one should add the so-called daisy (or ring) 
contribution 


501 


King {<p,T) = 


12t r 


(37) 


which represents the leading term of higher-order loop corrections that may play an important role during the 
EWPT dynamics. In V r i ng ((p,T) the sum is over the scalar and longitudinal gauge degrees of freedom, rh%((p,T) = 
+ Ili(T) are their thermal squared masses, and n^(T) are the thermal parts of the self energies, which are 
collected in appendix El To estimate V r i ng (ip,T), one performs the resummation of an infinite class of infrared- 
divergent multiloops diagrams, known as ring diagrams, that describes the dominant contribution of long distances 
and gives a significant contribution when (almost) massless states appear in the system. In our case, we will include 
this by following another approach. Rather than adding V T i ng (<p,T) to V e s(<p,T), we will replace in eq. (1351) the 
field-dependent masses of the scalar and longitudinal gauge degrees of freedom with their thermal masses fhi(<p,T). 
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*1 


FIG. 6 : The relative changes of the trilinear Higgs self-coupling (left), the Higgs mass (right) and the parameter (bottom) 
due to loop corrections versus the Higgs quartic self-coupling Ai. The Higgs mass rrih is fixed to 125.1 GeV. The blue points 
represent the SM values. 


In the criterion for a strong first-order phase transition, eq. (1341) . the critical temperature T c is the value at which 
the two minima of the effective potential are degenerate, 


d_ 

dtp 


V eS ((p,T c ) 


= 0, V e s(ip = v c ,T c ) = V eS (<p = 0,T C ). 


(38) 


In the SM, this leads to a Higgs mass below 42 GeV [5JLJ, since the ratio v c /T c is inversely proportional to the Higgs 


martic coupling Ai. The strength of the EWPT can be improved if new bosonic degrees of freedom are invoked [52- 
|55|, which is the case we are investigating. It is clear from eq. (l30l) that for large values of the couplings and/or masses 
of the extra scalars, the one-loop corrections to the Higgs mass could be significant, which allows Ai to be smaller 
and, therefore, fulfills the criterion in eq. (13411 without conflicting with the recent Higgs mass measurements 23[ 241. 
Here, the relevant couplings are those of the Higgs doublet to the charged scalars, A 3 a , and to the neutral ones, 
A 3 a + A 4 a , in the limit | A 5 1 <€C |A 3 a , 4 a|- The situation may be compared to those in similar setups 56-5^] where extra 
scalars can help bring about a strongly first-order EWPT by (a) relaxing the Higgs quartic coupling Ai to as small 
as 0 ( 1 O -4 ) and (b) enhancing the value of the effective potential at the wrong vacuum at the critical temperature 
without suppressing the ratio v c /T c , which relaxes the severe bound on the mass of the SM Higgs. 

The integral in eq. (1361) is often approximated by a high temperature expansion. However, in order to take into 
account the effect of all the (heavy and light) degrees of freedom, we will evaluate them numerically. 

With the same 5000 benchmark points used previously, in figure [3 we present v c /T c as a function of T c and of the 
Higgs quartic self-coupling. It is obvious that the criterion for a strongly first-order EWPT is easily satisfied for a 
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FIG. 7: Left panel: v c /T c versus T c , estimated with (red) and without (green) the daisy contribution. Right panel: v c /T c 
versus the Higgs quartic self-coupling Ai, estimated by considering the daisy contribution. 


large number of benchmarks. Moreover, we find that the daisy contribution to the effective potential tends to weaken 
the EWPT strength in this setup. One also notices that a strong EWPT can be obtained for different values of the 
Higgs quartic self-coupling Ai, as shown in the right panel of figure [3 even for values larger than the tree-level one, 
Ai. This leads us to conclude that the EWPT is always strongly first-order due the reason (b) mentioned above, 
where the extra heavy scalars’ existence makes the Higgs VEV slowly varying with respect to temperature and the 
wrong vacuum value, i.e. V e s(ip = 0, T), is evolving and increases with temperature. 

We remark that due to the absence of a CP-violating phase in the potential V an additional source of CP violation 
has to be included in the Lagrangian of the more complete theory for it to be realistic for baryogensis. One possibility 
is to introduce dimension-six operators which couple the inert scalars to the top-quark mass and are suppressed by a 
new-physics scale that can be well above one TeV, in analogy to a scenario of electroweak baryogenesis from a singlet 
scalar foil . 


VII. DISCUSSION & CONCLUSION 

According to the analysis carried out in previous sections, the extra scalars can have important effects on the Higgs 
phenomenology and the electroweak phase transition if these particles are relatively light and the couplings to the SM 
Higgs doublet are large (A 3 a for charged scalars and A 3 a + X^ a for neutral ones). Therefore, from the 5000 benchmark 
points used previously, we extract those that simultaneously satisfy (i) the constraint from the measurements on the 
Higgs decay mode h 77 , namely 0.9 < TZ 11 < 1.37, (ii) the electroweak precision tests, i.e. all the points inside the 
three ellipsoids in figured and (iii) the criterion v c /T c > 1 for strongly first-order EWPT. As mentioned in section llVl 
the Higgs decay channel into a pair of inert scalars is closed for all the viable benchmarks and hence its experimental 
bound is not relevant. Here, we divide the points fulfilling the conditions (i,ii,iii) into three sets according to the 
ellipsoid to which they belong on the (AS', AT) plane. The results are displayed in figure [H 

From the top panels in figure [5J one can see that the extra scalar masses do not exceed 900 GeV according to our 
parameter choices in eq. (fl3l) . The charged scalars could be light up to the LEP II bound (100 GeV), while the neutral 
scalars, which were supposed to be less constrained before, are now not allowed to be less than 120 GeV due to the 
electroweak precision tests in this model. From the bottom left panel, it is evident that the couplings of the Higgs 
doublet to the charged scalars, A3 a , and to the neutral ones, A3 a + A4 Q , could be both larger than 1 or smaller than 
0.5. They could vary also within the whole considered range [0:3], or they could be almost equal in absolute values 
(i.e. close to the dashed curve). The bottom right plot in this figure reveal that, while strongly first-order EWPT 
occurs for all of the viable benchmark points, only for some of them there is a positive correlation between the EWPT 
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FIG. 8: Top panels: the distribution of masses of the extra charged and neutral scalars under the assumptions described in the 
text. Bottom left panel: the strength distribution of the different quartic self-coupling of the Higgs doublet to the extra charged 
and neutral scalars; the dashed line represents the case where the couplings are equal in magnitude, i.e., Min(|A 3 a |, |A 3 a + A 4 a|) = 
Max(|A 3 a |, |A 3 a + A 4 Q |). Bottom right panel: the EWPT strength versus the relative enhancement of the Higgs trilinear self¬ 
coupling. The red, green, and blue points correspond to those inside the 99%, 95%, and 68% CL ellipsoids on the (AS, AT) 
plane in figure [4] 


strength and substantial enhancement of the Higgs trilinear self-coupling relative to the SM prediction as shown in 
ref. 


la¬ 


in conclusion, we have considered a scenario beyond the SM involving three scalar weak doublets and investigated 
a number of implications of the case where two of the doublets are inert and charged under a dark Abelian gauge 
symmetry. We looked at the effects of the new scalars on oblique electroweak parameters, the Higgs decay modes 
h -A 77 , 7 Z, and its trilinear coupling. We also examined how the inert scalars can induce strongly first-order EWPT. 
Taking into account various theoretical and experimental constraints, we demonstrated that the viable parameter 
space can all accommodate strongly first-order EWPT and contains regions in which the Higgs trilinear self-coupling 
is enhanced/reduced by up to 150% compared to its SM value. Future experiments with sufficient precision can test 
the new scalars’ effects that we have obtained on the Higgs decays h -A 77 , 7 Z and trilinear coupling. 
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Appendix A: Vacuum stability conditions 


We can rewrite the doublets $ and 771,2 and their products according to 


$ = /$, = 1 

&VaVl$ = Pa, 


vlmvlm = p‘ 


Va = e a rja, fllfla = 1 , f,e a > 0 , 

0 < Pa,p' < 1 - 


(Al) 


Assuming that A5 in eq. © is negligible compared to the other A’s, we can then express the part of V that is quartic 
in the doublets approximately as 


2„2 
2 


V4 — ^Ai/ 4 + ^A 2 ie 4 + iA 22 e2 + A31 / 2 e 2 + A3 2 / 2 e; 
+ Mif 2 ^ipi + Xi 2 f 2 e\p 2 + Xeefe^ + A-re^e^p' 

( P 


ei 4)X 


e? 


(A2) 


V 


where 


( Ai A31 + P1A41 A3 2 + p 2 X 4 2 \ 


A = 


A31+P1A41 A 2 i 

V ^32 + P 2 Xl 2 Xq + p' A7 


A6 + p'X 7 I . 

A22 / 


(A3) 


To ensure the stability of the vacuum, we need to derive relations among the A’s in V4, which dominates V at large 
fields, such that the minimum of V4 remains positive. This can be achieved using copositivity criteria [ 6 lJ , which in 
this case are applied to the minimum of A. Since A4 a ,7 can be positive, zero, or negative and 0 < p a ,p' < 1 , we have 


/ 


Amin — 


Al 


A31 + Min ( 0 , A41) A3 2 + Min ( 0 , Ai 2 ) 


A31 + Min ( 0 , A41) A 2 i Ag + Min ( 0 , A7) 

y A32 + Min ( 0 , A42) Ag + Min ( 0 , A7) A22 


(A4) 


From the criteria for strictly copositive 3 x 3 matrices 62 64 j then follow the conditions in eq. ©. 


Appendix B: Interaction terms for S a and V a 


The interaction terms of xi,2 in eqs. m and ( 1141 ) can be rewritten in terms of the real and imaginary components 
defined in eq. Q). Thus 


C D 


2 c w l 

+ f{ 


c 2 B (Vi - V 2 ^S 2 ) + s 2 e (v 1 pS 2 + V 2 7 f S^) 
csH+dPS 1 - iV 1) + c e H+dPS 2 + iV 2 ) 


+ ssH+ d 7 S 2 - iV 2) - ssH+d^Sx + iV 1) 

pi / fy l 

. _i_ W+Ml* ^ 


+ ^ (S? + V\ + + Vl) (j^- + W+»W~ 


W~- H.c.} 


h 


— [( c e^ 2 i + s eP 22 ~ m xi){Si + ^ 1 ) + ( c eP- 22 + s eP 2 i — m L) (^2 + ^ 2 ) 
+ s 28 (/i 2 i — p 22 ) {S\S 2 + ViV 2 )\ . 


(Bl) 
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Appendix C: Field-Dependent and Thermal Masses 

To estimate the Higgs effective potential, one needs the field-dependent squared masses m 2 {ip) of all the contributing 
particles. One also requires the first, second, and third derivatives of m 2 (tp) to determine the counterterm d/sf in 
eq. (l29l) . the one-loop correction to the Higgs mass, and the enhancement of the Higgs trilinear self-coupling. 

The field-dependent masses of the electroweak gauge bosons and top quark have their SM values. For the other 
particles, we have the thermal masses rhi = mi(tp,T) which are given by 

frig = Mi + + n$, rr? h = fi\ + |Ai ip 2 + n$, 

= g 2 a + 5 A 3 a ip 2 + n Va , m 2 Xl 2 = \ (Ci + C 2 T '/ r Rj, 

C’a = /X 2a + ^(A 3 a + A 4 a)y ) 2 + H^ a , R = (C\ — C 2 ) 2 + 4c 2 , C=||A 5 |<^ 2 , (Cl) 

and related to rrii(ip ) by m 2 (<p, T) = m 2 (ip) + n^, where H; = n*(T) denote the thermal parts of the self energies 
and n$ j7)a are listed below. Hence the Goldstone bosons (Q) and the Higgs boson also have the same field-dependent 
masses as their respective counterparts in the SM. We note that the inert CP-even and CP-odd neutral scalars mix, 
leading to equal-mass eigenstates, according to eq. ©. 

It is simple to get the first, second, and third derivatives of m 2 (<p) from eq. m- For completeness, here we supply 
them explicitly: 


rrig{tp) = \i<p, rh 2 h (ip) = 3Aivj, rh 2 H ±(ip) = \ 3a tp, 

-Op)-£(4 + 4^). 

Ca = (^30, + A 4 a) ^5 ^ = 2 

R = 2 (Ci - C 2 ) (Ci - C 2 ) + 8 cc, 


(C2) 


rh 2 g (ip) = \i, m 2 h (ip) = 3Ai, fh 2 H ±{ip) = \ 3a , 

2 M = \(ci + C* 


Xl,2 ' 


2 VR Ay/W I 


Ca — A3 a + A 4 q, C — A | A 5 1 , 


R = 2 


(Ci - C 2 ) 2 + 2 (C-, - C 2 ) (Ci - C 2 ) + 8 c 2 + 8 cc, 


(C3) 


™ gM = ™ 2 h(<p) = "4±(^) = Ca= c = 0, 

...2 , s 1 ( ... 3 RR 3 R 3 \ 

m xi,M = "F - 7 ^ I R- 


4 VR 2 R 4P 2 J ’ 

R = 6 (Ci - C 2 ) (Ci - C 2 ) + 24cc. (C4) 

Finally, we write down the thermal parts of the pertinent self-energies. For the scalar and electroweak bosons 


9 3 

n$ = ( 6A1 + — g 2 + —g 2 ' + 3 y 2 + 4A31 + 2A41 + 4A3 2 + 2A4 2 


T 2 

12 ’ 


9 2 ^ 2 


n r)a — ( —g + — gy + 4A3a + 2A4 a + 6A 2a + 4 Ag + 2A7 -I- 


3 ~2 

12 ’ 


„ 17 9 o „ 11 9 9 

n w = -frg~T , n B = — g\-T , 
6 16 


(C5) 
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where yt denotes the top-quark Yukawa coupling and Q Va = QcVa is the charge of the inert doublet r] a under U(l)i> 
Numerically, since g^> is unknown, for definiteness we set g^> = gy ■ 


[1] G. Aad et al. [ATLAS Collaboration], Phys. Lett. B 716, 1 (2012) [arXiv: 1207.7214 [hep-ex]]. 

[2] S. Chatrchyan et al. [CMS Collaboration], Phys. Lett. B 716, 30 (2012) [arXiv: 1207.7235 [hep-ex]]. 

[3] K.A. Olive et al. [Particle Data Group Collaboration], Chin. Phys. C 38, 090001 (2014). 

[4] J.F. Gunion, H.E. Haber, G.L. Kane, and S. Dawson, The Higgs Hunter’s Guide (Westview Press, Colorado, 2000). 

[5] G.C. Branco, P.M. Ferreira, L. Lavoura, M.N. Rebelo, M. Sher, and J.P. Silva, Phys. Rept. 516, 1 (2012) [arXiv: 1106.0034 
[hep-ph]]. 

[6] I.P. Ivanov, V. Keus, and E. Vdovin, J. Phys. A 45, 215201 (2012) [arXiv:1112.1660 [math-ph]]. 

[7] A.C.B. Machado and V. Pleitez, arXiv: 1205.0995 [hep-ph]. 

[8] V. Keus, S.F. King, and S. Moretti, JHEP 1401, 052 (2014) [arXiv:1310.8253 [hep-ph]]. 

[9] E.C.F.S. Fortes, A.C.B. Machado, J. Montano, and V. Pleitez, arXiv: 1408.0780 [hep-ph]. 

[10] M. Maniatis and O. Nachtmann, arXiv:1408.6833 [hep-ph], 

[11] I.P. Ivanov and C.C. Nishi, JHEP 1501, 021 (2015) [arXiv: 1410.6139 [hep-ph]]. 

[12] S. Moretti and K. Yagyu, arXiv:1501.06544 [hep-ph]. 

[13] V. Keus, S.F. King, S. Moretti, and D. Sokolowska, JHEP 1411, 016 (2014) [arXiv:1407.7859 [hep-ph]]. 

[14] E. Ma, I. Picek, and B. Radovcic, Phys. Lett. B 726, 744 (2013) [arXiv: 1308.5313 [hep-ph]]. 

[15] S. Kanemura, Y. Okada, and E. Senaha, Phys. Lett. B 606, 361 (2005) [hep-ph/0411354], 

[16] A. Noble and M. Perelstein, Phys. Rev. D 78, 063518 (2008) [arXiv:0711.3018 [hep-ph]]. 

[17] C.W. Chiang and T. Yamada, arXiv:1404.5182 [hep-ph]. 

[18] K. Fuyuto and E. Senaha, Phys. Rev. D 90, no. 1, 015015 (2014) [arXiv: 1406.0433 [hep-ph]]. 

[19] D. Curtin, P. Meade, and C.T. Yu, JHEP 1411, 127 (2014) [arXiv: 1409.0005 [hep-ph]]. 

[20] S. Kanemura, T. Kasai, and Y. Okada, Phys. Lett. B 471, 182 (1999) [hep-ph/9903289]. 

[21] S.Y. Ho and J. Tandean, Phys. Rev. D 89, 114025 (2014) [arXiv: 1312.0931 [hep-ph]]. 

[22] E. Ma, Phys. Rev. D 73, 077301 (2006) [hep-ph/0601225]. 

[23] G. Aad et al. [ATLAS Collaboration], Phys. Rev. D 90, 052004 (2014) [arXiv:1406.3827 [hep-ex]]. 

[24] CMS Collaboration, Report No. CMS-PAS-HIG-14-009, July 2014. 

[25] D. Curtin, R. Essig, S. Gori, P. Jaiswal, A. Katz, T. Liu, Z. Liu and D. McKeen et al., Phys. Rev. D 90, 075004 (2014) 
[arXiv: 1312.4992 [hep-ph]]. 

[26] A. Falkowski, F. Riva, and A. Urbano, JHEP 1311, 111 (2013) [arXiv: 1303.1812 [hep-ph]]. 

[27] P.P. Giardino, K. Kannike, I. Masina, M. Raidal, and A. Strumia, JHEP 1405, 046 (2014) [arXiv: 1303.3570 [hep-ph]]. 

[28] J. Ellis and T. You, JHEP 1306, 103 (2013) [arXiv: 1303.3879 [hep-ph]]. 

[29] G. Belanger, B. Dumont, U. Ellwanger, J.F. Gunion, and S. Kraml, Phys. Rev. D 88, 075008 (2013) [arXiv: 1306.2941 
[hep-ph]]. 

[30] K. Cheung, J. S. Lee and P. Y. Tseng, Phys. Rev. D 90 (2014) 9, 095009 [arXiv: 1407.8236 [hep-ph]]. 

[31] G. Aad et al. [ATLAS Collaboration], Phys. Rev. D 90 (2014) 11, 112015 [arXiv: 1408.7084 [hep-ex]]. 

[32] https://twiki.cern.ch/twiki/bin/view/LHCPhysics/CERNYellowReportPageBR3. 

[33] S.Y. Ho and J. Tandean, Phys. Rev. D 87, 095015 (2013) [arXiv: 1303.5700 [hep-ph]]. 

[34] C.S. Chen, C.Q. Ceng, D. Huang, and L.H. Tsai, Phys. Rev. D 87, 075019 (2013) [arXiv:1301.4694 [hep-ph]]. 

[35] B. Swiezewska and M. Krawczyk, Phys. Rev. D 88, no. 3, 035019 (2013) [arXiv: 1212.4100 [hep-ph]]. 

[36] A.D. Banik and D. Majumdar, Eur. Phys. J. C 74, no. 11, 3142 (2014) [arXiv: 1404.5840 [hep-ph]]. 

[37] M.E. Peskin and T. Takeuchi, Phys. Rev. D 46, 381 (1992). 

[38] M. Baak et al. [Glitter Group Collaboration], Eur. Phys. J. C 74, no. 9, 3046 (2014) [arXiv: 1407.3792 [hep-ph]]. 

[39] R. Barbieri, L.J. Hall, and V.S. Rychkov, Phys. Rev. D 74, 015007 (2006) [hep-ph/0603188]. 

[40] H. Baer, et al., ’Physics at the International Linear Collider’, available at: http://lcsim.org/papers/DBDPhysics.pdf 

[41] A. Ahriche, A. Arhrib and S. Nasri, JHEP 1402 (2014) 042 [arXiv:1309.5615 [hep-ph]]. 

[42] A.B. Lahanas and D.V. Nanopoulos, Phys. Rept. 145, 1 (1987). 

[43] S.P. Martin, Phys. Rev. D 65, 116003 (2002) [hep-ph/0111209]. 

[44] V.A. Kuzmin, V.A. Rubakov, and M.E. Shaposhnikov, Phys. Lett. B 155, 36 (1985). 



17 


[45] M. Trodden, Rev. Mod. Phys. 71, 1463 (1999) [hep-ph/9803479]. 

[46] M.E. Shaposhnikov, Nucl. Phys. B 287, 757 (1987). 

[47] M. E. Shaposhnikov, Nucl. Phys. B 299, 797 (1988). 

[48] L. Dolan and R. Jackiw, Phys. Rev. D 9, 3320 (1974). 

[49] S. Weinberg, Phys. Rev. D 9, 3357 (1974). 

[50] M.E. Carrington, Phys. Rev. D 45, 2933 (1992). 

[51] A.I. Bochkarev and M.E. Shaposhnikov, Mod. Phys. Lett. A 2, 417 (1987). 

[52] A. Ahriche and S. Nasri, Phys. Rev. D 83, 045032 (2011) [arXiv: 1008.3106 [hep-ph]]. 

[53] T.A. Chowdhury, M. Nemevsek, G. Senjanovic, and Y. Zhang, JCAP 1202, 029 (2012) [arXiv: 1110.5334 [hep-ph]]. 

[54] D. Borah and J.M. Cline, Phys. Rev. D 86, 055001 (2012) [arXiv:1204.4722 [hep-ph]]. 

[55] G. Gil, P. Chankowski and M. Krawczyk, Phys. Lett. B 717, 396 (2012) [arXiv:1207.0084 [hep-ph]]. 

[56] A. Ahriche and S. Nasri, JCAP 1307, 035 (2013) [arXiv: 1304.2055]. 

[57] A. Ahriche, Phys. Rev. D 75, 083522 (2007) [hep-ph/0701192], 

[58] A. Ahriche and S. Nasri, Phys. Rev. D 85, 093007 (2012) [arXiv:1201.4614 [hep-ph]]. 

[59] A. Ahriche, K.L. McDonald and S. Nasri, in preparation. 

[60] J.M. Cline and K. Kainulainen, JCAP 1301, 012 (2013) [arXiv:1210.4196 [hep-ph]]. 

[61] K. Kannike, Eur. Phys. J. C 72, 2093 (2012) [arXiv: 1205.3781 [hep-ph]]. 

[62] K. Hadeler, Linear Algebra Appl. 49, 79 (1983). 

[63] G. Chang and T.W. Sederberg, Comput. Aided Geom. Des. 11(1), 113-116 (1994). 

[64] L. Ping and F.Y. Yu, Linear Algebra Appl. 194, 109 (1993). 

[65] J.I. Kapusta and C. Gale, Finite-temperature field theory: Principles and applications (Cambridge University Press, Cam¬ 
bridge, 2006). 


